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I. INTRODUCTION 


In a previous paper 


Q 


(hereafter referred to as paper I), the equations of motion of com¬ 
pact binary systems were computed at the 3.5 post-Newtonian (PN) order, corresponding 
to the formal XjcJ level in a post-Newtonian expansion with respect to the Newtonian accel¬ 
eration (where c is the speed of light). This conhrmed and extended previous calculations 
in Refs. j2|, y, y, |5| . The 3.5PN terms consist of gravitational radiation reaction terms, IPN 
order beyond the lowest order of radiation reaction, which occurs at 2.5PN order in the ac¬ 
celeration. To deal with the divergences of Poisson or Poisson-like integrals in higher orders 
(a well-known problem in the standard post-Newtonian formalism), paper I used a particu¬ 
lar solution of the Poisson equation dehned by means of a hnite part operation resorting to 
complex analytic continuation. 

The aim of the present paper is to justify fully the procedure that was employed in 
paper I to compute high-order post-Newtonian terms (in particular the radiation reaction 
terms therein), and that was based on a formal expansion of the retardation of the standard 
retarded integral, augmented by the hnite part operation. The proof will essentially consist 
of deriving a new expression for the post-Newtonian expansion in the near zone of a general 
isolated post-Newtonian source, which will constitute an interesting variant of the solution 
proposed in a preceding article (hereafter paper II). The latter solution is obtained 
within a particular framework combining the multipolar-post-Minkowskian formalism with 
a matching to a post-Newtonian source [ll la, Isf- 

More precisely, in paper II, a general expression is provided for the gravitational held 
generated by an isolated (slowly-moving) system in its near zone, developed up to any post- 
Newtonian order. This expression is free of the well-known divergences of Poisson-type 


integrals which had plagued t 
earlier works (see e.g. Refs. 


le post -Newtonian expansion at high-order approximations in 

17| for discussions). The method 





adopted in paper II to deal with these divergences has been to replace the Poisson integrals 
by a new type of solution to the hierarchy of post-Newtonian equations. It involves a 
“regularization” of the boundary of the Poisson integral at inhnity by means of the hnite 
part. Let us emphasize that the hnite part is only used as a mathematical trick to hnd 
in a convenient way a particular solution of the diherential equation we want to solve. To 
this particular solution, we add a homogeneous solution which is determined by asymptotic 
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matching to the external held. It is then possible to iterate formally the post-Newtonian 
expansion ad infinitum. 


The matching must be understood as a variant of the method ^f 
matching applicable in the exterior part of the source’s near zone 


asymp 


la 


19, 


otic-expansion 


m, 


211. As the 


exterior solution satishes a condition of no-incoming radiation, the solution of paper II incor¬ 
porates the effects of radiation reaction and, therefore, gives a correct physical description 
of the isolated system. 

The alternative form that we derive for the end result of paper II further clarihes its 
relationship with the usual method of computing the successive post-Newtonian approxi¬ 
mations, based on a naive expansion of the held equations, like for instance in Ref. |l^ . 
Essentially, we show that the present method leads to correct results up to the 3.5PN order; 
this is, however, no longer the case starting from the 4PN order ~ 1/c®, because of the non¬ 


linear contributions of gravitational wave tails 0 I22L |23[ . As we shall see, the expression 


of the post-Newtonian expansion derived below provides further insights on the structure of 
the gravitational radiation reaction terms and the contribution of tails in the near zone. 

Section El recalls some necessary results from paper II, and ohers complementary com¬ 
ments. In Section EH we derive the novel expression for the post-Newtonian expansion, 
which we subsequently discuss. Finally, in Section Cvl we check directly, by using our new 
form, that this solution is divergenceless, and thus satishes the harmonic coordinate Einstein 
held equations. The paper ends in Section O with a short conclusion. 


II. REMINDERS OF EARLIER WORK AND COMMENTS 

In the following, we shall describe the gravitational held in harmonic coordinates, assumed 
to exist globally and to be unique, in terms of the “Gothic” metric deviation from hat space- 
time, We shall thus have 


duhf^’' = 0, (2.1a) 

(2.1b) 

where and g are the inverse matrix and the determinant of the usual covariant metric 
ga/B respectively, and where = diag(—1,1,1,1) denotes the Minkowski hat metric in 
Cartesian coordinates. With this choice and notation, the “relaxed” Einstein held equations 
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simply read 




IdvrG 




( 2 , 2 ) 


with □ = r]°‘^dai 3 being the flat d’Alembertian and G the Newtonian constant. The pseudo 
stress-energy tensor accounting for the matter as well as the gravitational fields in 
harmonic coordinates, is given by 

^4 


= (-g) + 


IQttG 




(2.3) 


where is the stress-energy tensor of the extended localized (bounded) matter source 
with spatially compact support, and is the gravitational stress-energy distribution, 

which depends non-linearly on and its space-time derivatives. The pseudo-tensor is 
conserved in the usual sense: 

= 0. (2.4) 


Now, we are interested in the general expression of the post-Newtonian expansion valid 
in the source’s near zone, he. for the coordinate distance r significantly smaller than the 
gravitational wave length of the emitted radiation. The held equations (EIII)-(E3), supple¬ 
mented by a condition of no-incoming radiation, are solved iteratively as a formal expansion 
when c —>■ -|-cx). The post-Newtonian series of the held variable denoted by means of 
an overbar, = PN , takes the following “polylogarithmic” form [3| : 

ft'“'(x,i,c) = 5^<hAAj;-(x,«), (2.5) 


with post-Newtonian coefficients (p, g G N) depending on the field point x and coordinate 
time t. The pseudo-tensor, admits a post-Newtonian expansion of the same type, 
say with coefficients The post-Newtonian expansions and are connected via the 
field equations by 


dV = ^7“-. 


( 2 . 6 ) 


completed with the harmonic gauge condition 


dX’' = 0 . 


(2.7) 


In the present paper, we shall not consider expressions of the post-Newtonian coefficients 
at a given order, but rather work with whole post-Newtonian expansions and of 
the type (EH). Such expansion series will be manipulated by means of the rules of formal 
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series. In particular, we shall not investigate the eventual convergence properties, or, more 
generally, the mathematical nature of these series. The developments presented in paper II 
and the present one will thus constitute some formal generalization, to any post-Newtonian 
order, of the expressions of the near zone metric and the radiation reaction force. These 


computations at 2PN order for the 
and at 1.5PN relative order for the 


expressions have been explicitly validated by detail^ 
matching process as well as the multipole moments 
radiation reaction contributions j22, |23| . 

The strategy adopted in paper II was (i) to solve the wave equation (ESD in the form of a 
certain functional of the pseudo tensor, , and (ii) to verify that the gauge condition 

(ITTIl is satished as a consequence of the conservation of Eq. (El). Importantly, the 
post-Newtonian expansion of the gravitational held, was found to be a functional not 
only of the post-Newtonian expansion of the pseudo-tensor, , but also of its multipole 
expansion, denoted below by which is valid in the source’s exterior, and in particular 

in the asymptotic regions at inhnity. This can be understood by the fact that the post- 
Newtonian solution (comprising the radiation reaction terms) depends on the boundary 
conditions imposed at inhnity, especially the no-incoming radiation condition at past null 
inhnity, which must describe a matter system isolated from the rest of the universe. 

In the formalism of paper II, the no-incoming radiation condition is trivially satished, 
because it is assumed that the matter source is stationary before a hnite date in the remote 
past, say, for t < —T, with —T denoting a hxed instant in the past. More precisely, we 
suppose that both and are constant (with respect to time) outside the causal future 
development of the source when it starts emitting radiation, i.e. for instance 


r^^(x,f) = rr(x) when t-'-^<-T, 


( 2 , 8 ) 


where the function rg^^(x) does not depend on t, and refers to the “stationary” value taken 
by the pseudo stress-energy tensor at early times. In the present paper, we shall continue 
to assume that Eq. (EHD and the same equation for hold. It will be convenient to 

decompose the pseudo tensor into the sum of its stationary value, dehned by Eq. (ESD, and 
a “non-stationary” part, denoted by r^J^(x, t), that vanishes when t < —T. Hence, we pose, 
at any time t, 

r'^^(x,f) = rr(x) + rr(x,f). (2.9) 

Such splitting is always possible thanks to the hypothesis of past stationarity. In this 
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formalism, the derivation of any physical result is to be done under the assumption (EHD, 
but at the end of the calculation, once the result is in hand, T can be sent to +cxd. 

The main outcome of paper II is that the post-Newtonian solution of the wave equation 
(Eni), and satisfying (Eza), together with correct boundary conditions at inhnity, may be 
written as the sum of a particular (inhomogeneous) solution plus a specihc homogeneous 
solution; both of which are now described. 

Let us hrst consider the particular solution. It is given by the so-called operator of 
instantaneous potentials, denoted by which is dehned by a straightforward iteration 
of the inverse Poisson operator A“^. The instantaneous operator is endowed with 
a regularization, dealing with the boundary of the integrals at inhnity, i.e. for |x'| — 
-|-cx). This regularization is dehned from a particular finite part procedure, indicated by 
an operator FPb =05 with analytic continuation in the parameter B E C down to i? = 0. 
Essentially, it consists in multiplying the integrand of the Poisson integral by a factor |x'|^, 
and taking the hnite part of the Laurent expansion of the (analytic continuation of the) 
integral when B —y 0. The analytic continuation is always possible at a given post-Newtonian 
order for the particular class of functions considered in paper II as well as in the present 
paper (c/. Appendix B in paper II). Thanks to this particular regularization, ^ we can select 
a solution of the iterated Poisson equations, that is diherent from the standard Poisson-type 
integral and avoids the appearance of divergences at the boundary at inhnity (see paper II 
for the details). The crucial point is that the latter construction is a particular solution of 
these equations. We thus have 

m\Tn = Y.[73t) ( 2 . 

fc =0 ^ ^ 



where the fc-th iterated Poisson integral reads 

"d^x' lx 


J\2k-l 


A-fc-i [rH (x,t) = 




( 2 . 11 ) 


—dvr (2/c)! 

The overtilde on a quantity means here (following the notation of paper II) regularization 
by use of the hnite part operation FPb=o: 


A-^-i [/(x)] ^ FP A-^-i [|x/ro|^/(x)] , (2.12a) 

£ 5=0 


^ Note that although it is sometimes convenient to call it as such, the scheme is strictly speaking not a 
regularization since there is no assumption that some infinite quantity should be formally discarded. 
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J d^x.' g{x.') = FT J rfV |x7ro|^ c/(x'). (2.12b) 

Following earlier articles (Refs. and paper II), we introduce an independent length 

scale tq inside the regularization factor in order to adimensionalize it. Anyway, we know 
that the constant tq disappears from the physical metric as soon as the metric is functionally 
related to the local matter distribution. As shown in paper II, Eq. dnni) indeed constitutes 
a solution of the d’Alembertian equation (EH): 


□ <! J-i Irn !> = 


(2.13) 


We next add a homogeneous solution in such a way that the post-Newtonian metric 
matches the external one, which is globally defined outside the source in the form of a mul¬ 
tipolar expansion. Clearly, such a homogeneous solution must be regular inside the matter 
source, and in particular for r —> 0, assuming that r = 0 corresponds to the spatial origin 
with respect to which the multipolar expansion is dehned. Therefore, the homogeneous so¬ 
lution must be made up of a sum of retarded (multipolar) waves minus the corresponding 
advanced waves; it is well known that such anti-symmetric waves are regular in a neigh¬ 
borhood of r = 0, and are associated with radiation reaction effects. ^ In the Fourier 
representation, their contribution takes the form of a superposition of plane monochromatic 
waves (see below). However, following paper II, we prefer a parametrization based on certain 
multipole-moment functions in the physical domain. 

In paper II, we showed that, when the particular solution is taken to be that of the instan¬ 
taneous operator, defined above by Eq. (ITTUl) . the multipole-moment functions parametriz¬ 
ing the homogeneous solution, say A^(t), split in fact into a sum of two terms, ^ 


yif (t) = Tf (t) -47r3^f (t). 


(2.14) 


The hrst of these functions, 3^j^’^{t), turns out to be equal to the actual source multipole 


^ Here, the qualifier “anti-symmetric” refers to the formal structure of the waves, i.e. retarded minus 
advanced, but not to their real behavior under a time reversal, which can be very complicated in high 
post-Newtonian approximations. 

^ Here, L = ii ■ ■ ■ ie stands for a multi-index composed of £ spatial indices {£ being often referred to as the 
multipolar order). For the second term in Eq. 12.1411 . we employ a slightly different notation than in paper 
H, where it is given by (paper H) = 
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moment of the isolated system, as seen from the exterior region [9|. It is given by 


ion 


-+1 


x'l / dz 6i{z) {x.',t — z\x.'\/c), (2.15) 




with x'j^ = STF ■ ■ -x'J denoting the symmetric-trace-free (STF) part of the product of 
i spatial vectors; the function d£{z), whose integral is normalized to unity, is dehned by 

( 2 £ + 1 )!! 


SAz) = 




2^+H\ 


'•+1 


dzSAz) = 1 . 


(2.16a) 

(2.16b) 


'-1 


We emphasize that expression (ITT31) . together with (ITTBll . is physically valid only for post- 
Newtonian sources. As such, it must be considered in a perturbative sense and given as some 
formal post-Newtonian expansion. The explicit post-Newtonian series reads (see paper II) 


+ CO 

?r(«) = E 

j=i) 


IXx'A-Hx'OT>“'(,x',t), 


(2.17) 


where we have introduced the convenient notation 

^ 2ij\{2i + 2j+ 1)\\ 




Li 


(2.18) 


which is justihed by the fact that A-^ [A (x'^)] = 

As can be seen from Eqs. (ITT3D or (ITT7I) . the function is dehned as a functional of 
the post-Newtonian expansion of the pseudo-tensor valid in the source’s near zone. The 

second piece of the multipolar function (EH, namely , directly issues from the matching 
process to the external held, and arises purely from non-linearities. Thus, it depends on the 
expression of the pseudo-tensor in the exterior of the source (including the region located at 
inhnity), where it is given in the form of a multipolar series Al(r^^), with M representing 
the multipolar expansion. Actually, in the present formalism, and Al(r^'^) as well, 

are identified with the so-called multipolar-post-Minkowskian expansion of the external held, 
in the sense of Ref. ^ The function may be written in a compact way as (see 221 


^ The post-Minkowskian expansion (or non-linearity expansion when G —*■ 0) plays a crucial role in the 
construction of the multipolar-post-Minkowskian solution in Ref. |3, but shall not be considered for 
the present purpose. Indeed, we shall assume here that we are always entitled to sum up formally the 
post-Minkowskian series and to treat the formal limit as if it were “exact”. 
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and paper II) 


= J dz'ye{z)M{T^){x',t-z\x'\/c), (2.19) 

where ji^z) = —2Se{z) is such that 

+ CXD 

dz-fi{z) = l. (2.20) 

Remark that the above normalization of the integral of jilz) is a consequence of that of 
6i{z) [Eq. (j2.16bjl ]. together with the fact that ” dzji{z) = 0 by analytic continuation in 
the variable £ G C. 

Note that, as indicated in Eq. the function involves the non-stationary part 

of the pseudo-tensor, given by Eqs. (EHD-Q; hence the integral over z is well-dehned, 
as it ranges in fact on a hnite domain [1, Zmax], where Zmax -l- 1 = c (t -|- T)/|x'|. This means 
that “31^ depends on the initial instant of emission, —T, but since it appears only through 
an anti-symmetric type wave, in Eq. (HIl) below, the dependence on —T hnally cancels out. 

One can show that an alternative form for is 

x'l/c) I , (2.21) 




' dz‘:i 


-dvr 


1 (- 2 ^- 1 ) 



where d'j^ is the STF multi-derivative operator associated with the integration point x', and 
the superscript (—2£ — 1) over denotes the (2£ -f l)-th time anti-derivative that 

vanishes in the remote past, say, when t < —T. The STF spatial multi-derivative admits 
the following explicit expression: 

'f{t — r/c)\ . {^-\-i)\ fd-3(^t — r/c) 


dr 




^^2T!(f' —i)! c3 * 

i=0 ^ ' 


( 2 . 22 ) 


where r = |x|, hi = = STF {di^ ■ ■ -di^)] the superscript {^ — i) means the (£ —z)-th 

time derivative. 

With being the sum of these functions as given by Eq. (12.141) . the post-Newtonian 
expansion in the form obtained using the method of matched asymptotic expansions in paper 
II reads: 

4G ^ (-1)^ A \ A^]^{t - r/c) - A^^{t + r/c) 

^ ^ A ^ d. M 2r 

1=0 ^ 

The explicit expression of the multipolar waves can be written from Eq. (j2.22j) (with c —*• —c 
in the case of the advanced wave), but we shall give some alternative formulas in Eqs. (ED) 
and below. 
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III. NEW EXPRESSION OF THE POST-NEWTONIAN EXPANSION 


The idea for determining the new form of the post-Newtonian expansion is based on the 
remark that, in Eq. the part of the multipolar anti-symmetric wave that is associated 

with the function corresponds in fact to the “odd” parity part (he., the part composed of 
odd powers of 1/c) of the retardation expansion in the standard retarded integral. Since the 
operator of the instantaneous potentials (ITTU1) - (ITTT1) precisely equates the corresponding 
“even” part, we can reconstitute the complete retarded integral of the post-Newtonian source 
term in a nice way. 

To begin with, we notice that the post-Newtonian expansion of the multipolar anti¬ 
symmetric wave reads 


dr 


f Tf (t - r/c) - Tf (t + r/c) ) _ 
\ 2r 


+ 00 


(2^+1)” to 


(xl) 


d 


— 
cdt J 


\ 2p+2l+l 


(3.1) 


with notation (ITTHl) . Next, we recall that the function itself is given by the expression 

()2.15|) or alternatively by the post-Newtonian series (ITT7I) . After inserting Eq. (ITT7I) into 
(ED), we are led to a rather complicated double expansion. Furthermore, we have to sum 
up the latter result over all multipolarities £ = 0, • • • , -fcx). At this point, let us change in 
an appropriate way the summation indices (posing in particular k = j + p + €). We are then 
able to recognize a double summation which turns out to be nothing but the expression of 
the quantity |x — x'p^, when developed in terms of STF harmonics. It can be written in a 
particularly economic form according to the following formula, in which we resort again to 
notation (ITTSD : 


J\2k 


{2k + l)\ 


E 


-1)' 


k-t 


£l (2£+l)!! ^ 

£=0 ^ ' p =0 


(xl) A- 


(X, 


(3.2) 


This formula can be viewed as a finite expansion in terms of Legendre polynomials, since 
the product A~p (xl) (L^) (where p, g G N) is proportional to xlx'^^, and the latter is 
known to be in turn proportional to the £-th Legendre polynomial of the scalar products of 
the unit vectors n = x/r and n' = x'/r', namely Pe (n ■ n'). Thanks to Eq. ()3.2j) . we arrive 
at the interesting relation 


-l-CXD 


(-1)' a f jr(t- r/c)-jf(i + r/c) l _ 




e=o 




V 


2r 


^ \cdt J 

k=0 ^ 


d^x' 


,/\2k 


{2k + l)\ 


(3.3) 
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which shows, as announced before, that the anti-symmetric wave associated with is 

equal to the odd part [i.e. corresponding to n = 2/c -|- 1, fc G N, in Eq. (Id .411 below] of the 
expansion of the retarded integral as regularized by means of the hnite part FP. 

In order to gather terms (ITTUl) and o, which are both present in Eq. we need 

to consider the new object that is made of the formal expansion of the standard retarded 
(R) integral but where each of the terms is regularized by means of the FP procedure 
(indicated by an overtilde) to deal with the divergences appearing on the boundary at 
inhnity: 

5 ? im (x.<). £ |x-xr-‘T-(x'.t). (3.4) 

This object should be carefully distinguished from, say, involving the retarded 

integral without any kind of regularization and with the retardations being kept un¬ 
expanded inside. Eq. (Q is also a priori different from the more “global” solution 
/(WtiG), in which the pseudo-tensor has not been expanded in post- 
Newtonian fashion. As we shall see below, Eq. (D differs also from the post-Newtonian 
solution itself, i.e. h^^/(167rG). 

Let us emphasize that Eq. (D constitutes merely the definition of a formal Taylor-type 
post-Newtonian expansion, each term of which is built from the (also formal) post-Newtonian 
expansion of the pseudo-tensor. Such a dehnition is of interest as it corresponds to what one 
would intuitively think is the “natural” way of performing the post-Newtonian iteration, i.e. 
by Taylor expanding the retardations. Moreover, each of the terms of the series dH are 
mathematically well dehned thanks to the FP procesS;^nd can therefore be implemented in 
practical computations like in the one reported in Ref. [l|. The important point, and the only 
one which matters for our purpose, is that the dehnition dH solves, in a post-Newtonian 
sense, the required wave equation, namely 



and constitutes therefore a legitimate particular solution for determining the general form 
of the post-Newtonian expansion as the sum of this quantity and some homogeneous 
solution which is necessarily of the anti-symmetric type. 

Combining Eq. dSH with our dehnition of the operator of the instantaneous potentials 
fl^ - (ITTTl) . which corresponds to the even part of the retardation expansion [i.e. n = 2k, 
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A; e M, in Eq. (13.411 ] . we find 


E? [T-i. FI [± J + 

£=0 ■ 


(3.6) 


This formula is the basis of our writing of the new form of the post-Newtonian expansion. 
Indeed, by comparing Eqs. and (uni), we obtain the following expression: 


IQttG 
h = —:— 


□r' 


T 


■pv 


+ 00 

£=0 


- 1)^ 5 \ {t - r / c) -{t + r / c) 
£l 2r 


(3.7) 


which relates the post-Newtonian expansion of = (WttG/c'^) to the developed 

operator acting on the post-Newtonian expansion of the source term, Eq. (I3.7jl . 

together with the conservation law for the pseudo-tensor [Eq. (12.411 ]. constitute our final 
solution for the present problem. 

Let us insist that the formula (jnii should be understood, for practical calculations, in 
an iterative post-Newtonian way. At the lowest order, the source has compact support 
and there is no need of the finite part procedure. Higher-order equations are obtained by 
inserting previous iterations into the source term [or Ad(r^^)] truncated at this order; the 
source will be of non-compact support because of the gravitational field contribution. One 
then multiplies it by the regularization factor |x/ro|^, integrates and performs the Laurent 
expansion when H 0 to get the finite parts needed in the two terms of Eq. (HIl). Some 
examples of such finite parts of integrals can be found in paper I [see e.g. Eqs. (4.15-16) 
there]. 

Equation (HID is more interesting and probably more “fruitful” in applications than the 
old form ()2.23|1 . Indeed, recall that the anti-symmetric waves in Eq. ()2.23j) . parametrized by 
AM^' _ are linked with radiation reaction effects. It was noticed in paper II that 

the part associated with starts at the usual dominant order for the radiation reaction 
force, i.e. the 2.5PN order (where it yields the standard expression for the “Newtonian”, 
quadrupolar radiation reaction). It includes also the next term of order 3.5PN corresponding 
to the mass octupole and current quadrupole radiation reaction, and probably also the sub¬ 
sequent contributions of higher multipole moments that parametrize the radiation reaction 
at the linearized level. ^ On the other hand, the part linked with the functions has been 


5 


By multipole moments here, we mean the ^‘’source’’’ moments defined in Ref. [9j, which couple together 
non-linearly in the expressions of the radiation reaction force and the “radiatife” moments seen at infinity. 
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shown in Ref. 221 to be due to the non-linear contribution of the gravitational wave tails 

n 

in the radiation reaction force. It was determined beforehand in Ref. [20|| and shown there 
to appear at the 4PN order, i.e. 1.5PN order relative to the dominant radiation reaction 
effect. 

Consequently, when computing the radiation reaction force limited at the 3.5PN order, 
the part associated with can be neglected, and we can keep only the part given by the 
hrst term of Eq. (EH); this corresponds indeed to the intuitive way of performing the post- 
Newtonian iteration (bypassing the tail effects), which was advocated in old investigations 
of the post-Newtonian expansion for isolated fluids, like those of Ref. [l^. Thus, Eq. ()3.7|) 
rigorously proves that the 3.5PN radiation reaction must be computed by considering the 
odd-parity part of the retardations in the usual retarded integral, provided that the source 
terms be multiplied by the regularization factor |x'|^ and the integral computed fo 
the FP prescription when R ^ 0. This is exactly what has been done in the paper 
notably Section IV there) devoted to the computation of the 3.5PN radiation reaction in the 
case of compact binary systems. However, for higher-order computations, at the 4PN level 
and beyond, it is definitely necessary to take into account the contribution of the functions 
in Eq. (13.7j) which is induced by wave tails. 

For some applications, it may be useful to rephrase our result with the help of some 
more systematic notation. Notably, we may use an appropriate operator, denoted 
henceforth, which maps the source term of the Einstein field equations, namely the pseudo 
tensor onto the post-Newtonian solution given by Eq. (HID- Let us first recall that = 
{IQtiG/c^) which allows us to suppress the overall factor IbvrG/c^ in Eq. ()3.7j) and 

reexpress the latter equation in the form of 


lowing 
l| (see 




(3.8) 


where the specihc homogeneous, anti-symmetric wave type solution is defined by 


+ 00 






(-1)^ b f - ‘Kit + r/c) 

-2;- 


(3.9) 


Eq. (13.8j) suggests introducing some convenient functional operators acting on the compo¬ 


nents of For the first term of Eq. (EHl), we introduce the operator whose action 
on is defined by 

(3.10) 
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With this notation, we include into the dehnition of both the FP process, indicated by 
the overtilde [see Eq. (ITTl) ]. and the operation of taking the post-Newtonian expansion, 
= PN[r^'^]. For the second term of Eq. ()d.8|l . we pose and hence FC 

represents the particular functional that is obtained when the functions in Eq. (El are 
replaced by their explicit expressions in terms of as given by Eqs. or (I2.21|) . Note 

that and 'K^^' depend in fact on the multipole expansion of the pseudo tensor, 

rather than on itself. As a result, the operation of taking the multipole expansion {M.) is 
included into the above dehnition of the functional Naturally, the latter functional is 

also a function of the held point x and of (coordinate) time t. Thus, our notation means more 
precisely h{[r^'^](x, t) = h{^'^(x, f). In a similar way, it is convenient to make use of what we 
call Oil, namely a functional of and a function of time t dehned by 
We can then rewrite Eq. (EH as 

□-1 [rM^] = [r^^] , (3.11) 

where the operator is naturally split into the two terms appearing in Eq. dH) as 

p-^ = n^ + F{. (3.12) 

We shall use this operator in the next Section devoted to the harmonicity conditions. 

Note that the anti-symmetric multipolar waves such as !K^^, given in Eq. dSSD, can be 
represented as a Fourier superposition of plane waves. We give here the relevant formulas 
without proof. The Fourier representation is arguably simpler and perhaps more physically 
appealing, but the alternative STF multipolar anti-symmetric formulation dSH, which we 
mostly use in the present work, displays a structure that is often more useful in applications. 
Both forms are linked by 

h{^"(x, t) = j d^k [A^^t(k) + A^^{k) (3.13) 

where the Fourier amplitudes ^((((^(k) and A([['(k) read, in terms of the multipolar functions 
fk^^(f), as 

(3.14a) 

* e=o 

^ mtK-ck), (3.14b) 

* e=o 
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with dt (t) representing the one-dimensional Fourier transform 

of Here, we pose k = |k|, = /cjj ■ ■ ■ ki^ and the brackets surrounding indices denote 

the STF projection. Conversely, the function is entirely determined by means of the 

Fourier amplitudes as® 

Sf(i) = (2<+1)!! j #k + A^{k) . (3.16) 

IV. HARMONICITY CONDITION 

In this Section, we prove that our post-Newtonian solution (EUl) is divergenceless as a 
consequence of the harmonicity condition, and thus constitutes an actual solution of the 
Einstein held equations. More precisely, we check that the conservation of the pseudo¬ 
tensor, Eq. (El, implies the satisfaction of the harmonicity conditions (EH- This proof 
has already been done in paper II, but we present here the main steps of an alternative 
demonstration, using some of the tools employed there, and based directly on the new form 
dSU of the post-Newtonian expansion. 

In order to verify the harmonicity condition, it is convenient to establish a sufficient 
property involving the operator T>~^ introduced in Eq. dSH- In fact, the solution written 
in the form (IXTTl) will be divergenceless as a consequence of = 0, if (mathematically 

speaking) we are allowed to commute the partial derivative with the 'D~^ operator, he. 
duV~^ = V~^ = 0. Let us define the following commutator: \du, V~^'\ = 

du V~^ — V~^ dy. Thus, it is sufficient to prove that 

V~^] = 0. (4.1) 

The commutator corresponding to the partial time derivative is seen to be zero, 
\d/dt, = 0, because the operations of taking the spatial integral and performing spatial 
differentiations which are involved in 'D~^ are obviously transparent to the time derivative. 

All what remains for us to prove is that = 0, for the partial space derivative 

® This relation is easily checked to be true by virtue of the formula 

y^tj nyt) = i^e+m-iY+ydL:Kn{o,t), ( 3 . 15 ) 

where (x, t) is given by (13.911 . 
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di = d/dx^. Therefore, we shall now show that 


[a., M] = -[di, ns*]. 


(4.2) 


The hrst stage of our proof consists in evaluating by decomposing the multi¬ 
derivative acting on the anti-symmetric wave {t — r/c) —{t + r/c)] into 

symmetric-trace-free (STF) pieces using the technical formula 


didh = diL + 


2i+l 




(4.3) 


where the brackets {■■■) denote the STF operation for the indices they enclose. Each 
multipole order ^ in the right-hand side of Eq. (O is thus given by a sum of two terms. 
After changing appropriately the summation indices and bearing in mind that satishes 
the homogeneous wave equation, so that = {c~^d/dt^'K^, we arrive at 


Y f 1 


£=0 


£l 


-h r/c) 


2£ + 3 


) 2r 
\ 2 

j + r/c) — (^idem for the retarded 


wavej 


(4.4) 


In a second stage, we shall work on the quantity ‘K[diT^'^] in order to obtain an alternative 
expression involving rather than The idea is to integrate by parts the 

integral on x' dehning as ^ 


r 

-dlLWn = FF -- 
B=o / -dvr 


b+CXD 




dz'yeiz) M{diT^^){x.'Y - ^|x'|/c), (4.5) 


in a way that the source term does not contain any more derivatives. To begin with, we 
observe that Ai{diT^^) is equal to a total space derivative, plus a term that turns out to 
involve a time derivative: 

VfWT-rXx'.t-zlx'l/c) = Si[A4K,‘')(x'.t-z|x'|/c)] + zn[ A4(t"‘')(x', i - z|x'|/c). 

(4.6) 

The hrst contribution may be integrated by parts, which has the effect of transferring the 
operator di onto the product |x'|^a;)^. Its explicit action on the latter quantity reads 


^idx'l £)^) = [Bn[h'L + {2£ + l){n[h'^ - n[^)] |x' 




(4.7) 


Henceforth we pose rg = 1. 
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As = 0 in the remote past, the surface term does not contribute. By applying 
the chain rule, we rewrite the second contribution in the right-hand side of Eq. (gSl) as 
— (za:'/|x'p)(c}/5z)Af(T^g'^)(x',t — z|x'|/c). In particular, the contribution to the corre¬ 
sponding integral on 2 ; may be integrated by parts in turn, and hence we are led to consider 
the first derivative (d/dz) (^ 7 ^( 2 ;)), which verihes {d/dz) = {2l-\-l){j(^{z) — 7 £_i( 2 ;)). 

As a result, we get 

/* ^+00 

*t|Sir'“](i) = |P y y dz[-Bn'Alt(z) 

+ (2f + 1)[(7H«) -7 »-i(z))>'4 + 7(-i(z)(>' 4 - nX)]}v<(rJ‘")(x',( - z|x'|/c). (4.8) 

The crucial point is that only a pole cx 1/5, generated by the radial integration, may give 
rise to a non-vanishing contribution for a source term proportional to B in the limit 5 — 0. 
Therefore, we may restrict the integration domain over the hrst term in the curly brackets 
above to a ball of radius R centered at the origin without affecting its value. An immediate 
consequence is that can be replaced there by its post-Newtonian or near zone value 

M-ijns)- Moreover, we have A4{tqs) = Af(r(^/) by the matching between the inner and 
exterior helds (see paper II). 

The second term entering the brackets of Eq. (USD, in view of the relation 
{d/dz)^{ji+i{z)) = {2i + 1 )( 2 £ + 3 )( 7 £_i( 2 ;) - 7 ^( 2 ;)), is nothing but the second derivative 
of 7^+1 (; 2 ) up to a constant factor. This derivative is handled by means of two successive 
integrations by parts, and the resulting derivatives that bear on are transformed 

like ^ ^ 

(£) (J/) l/c). (4.9) 

Taking also into account the fact that h'^ — n'h/ = —-1- 1), this yields 

r j 3 „/ r+00 

= y I dz^e{z)M{T>:^){^',t-z\^yc), (4.10) 

which may be used to eliminate from the anti-symmetric wave TC[c?jr^^]. 

At this level, we note that at any hnite post-Newtonian order, the structure of Af(r((/) 
is such (see paper II) that it yields a radial integral in the right-hand side of Eq. ()4.1()j) . 
that is decomposed into a sum of terms of the type dr'r'^^^ with g G Z. Now, a 
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useful lemma, already employed in Refs. 


and paper II, states that = 0 


by analytic continuation in B. This implies that the right-hand side of Eq. dnn will give 
zero if we make the substitution which means that can be replaced there by 

the complete pseudo tensor, -f Our desired commutator, i.e. the left-hand 

side of Eq. (14.211 . follows immediately from Eqs. (El) and ()4.1U|1 . We hnd 


[9..«]Ki = f 

e=o 


dr 


S^[rn(t-r/c)-Sl[r'^1(f + r/c) 


2r 


(4.11) 


where 


§*, Wn = FF B 


B=0 


—47r 


^H-oo 




dz^i{z) t — z\'x!\lc). 


(4.12) 


In the last stage of our proof, we shall focus on the other commutator in the right-hand- 
side of Eq. (US, which is readily obtained by integration by parts as 


+ CXD 




fc =0 


FPR 

B=0 


—dvr 




X X 


(4.13) 


We shall put it under a form where it is immediate to see that it is the opposite of TC], 
which will complete our demonstration. Due to the presence of a R factor, the only con¬ 
tributing terms are those that generate poles in the variable B due to the divergent behavior 
of when |x'| —>■ -|-cx). Thus, the integral depends only on the exterior domain, which can 
be chosen to lie outside a sphere of arbitrary radius R. If R is large enough, we are allowed 
to replace |x — and by their multipole expansions, so that 


FP B 

B=0 




'\B-2i\ 


—dTT 

FPR 

B=0 


xjx 




\x'\>R 


—dvr 


s'|'’-"i'M(|x-x'|‘-‘)M(r'“')(x', i). 


(4.14) 


By virtue of the lemma mentioned above, we can change the integration over the domain 
|x'| > R into an integration over the complementary domain |x'| < i? if we also change the 
sign. On the other hand, the expression of AI(|x — x'|*) for s G N is known explicitly. It is 
derived from the STF decomposition of the factors x'j^ (see e.g. Eq. (2.2) in j2^ or (A20a) 


m 


0) 


that enter the term of order i in the series expansion of |x — x'|^ valid for |x'| < |x| 
We recover the same formula as Eq. (n but for k = s/2, namely 

+ 00 / ^ +00 


X — X 


(s + 1)! 


E 


(-1)^ 


i\ (2 £+l)!!^ 
e=o ^ ’ p=o 


^ A-P { xl ) a 


-h+i+v 


A), 


(4.15) 
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provided that we define, for any integer s, 


-1- li" 

which generalizes the notation (ITTHD . We can transform (irail . with the help of the Euler 
function B(p, q) = 2 dz (z^ — 1)'^“^, to 

(—l)^+^( 2 f' + 1 )!! 


A 2 (xl) = 


■B(-e-(s + l)/2, e + i). 


(4.17) 


s! 2^+W. 

The appearance of the Euler function is of particular interest for our purpose, as it can simply 
be represented by an integral involving j^lz). Indeed, we have B(—£ — (s + l)/2, £ + l) = 
2 dzz^{z‘^ — 1)^, which immediately gives 


A-Hxx) = ^-^ dz^,{z)z\ 


s\ 


(4.18) 


'1 


Inserting the value of Al(|x—x'|*^“^) deduced from Eqs. ()4.15|1 - ()4.18|1 into the spatial integral 
with the radial variable |x'| ranging from 0 to i? (taking into account the change of sign 
when going from the outer to the inner domain), we obtain 

_ ^'^k+l _ 


[ft, DJ'] IrT = - ^ 


fc,£,p>0 


( 2 £+ !)!!£! 


A-P ( xl ) 


X I 4- 1 FPB 

cdt) B =0 


dz'yd 

'|<ij-4vr {k-2^-2p-l)\ 4 


^+oo 


dz-fi{z) , t). 

(4.19) 


Note that only the terms for which k > 2i+2p+l contribute to the sum over k since otherwise 
l/(fc — 2£ — 2p — 1)! = 0. The divergences near |x'| = 0 can be cured by means of the finite 
part regularization FP. We next permute the summation over k and the integration on 2 ;. We 
easily recognize, from the summation over k, the Taylor expansion of Al(rP^)(x', t — 2 ;|x'|/c) 
when c —> +cx). Finally, we obtain the relation 

2£+2p+l p r +00 

EPS / — —\x!\^~‘^x[x'i^ / dz^i{z) — z\y!\/c). (4.20) 

^=0 J|x'|<R -4vr 

Thanks to the formula (EH), which is valid for any function, we identify the right-hand 
side of with the opposite of the commutator found in Eqs. when it is 

developed as c —>• -|-cx). We can therefore conclude that 



v-^] = Dr'] = 0 , 


(4.21) 
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which, as we discussed above, achieves the proof that satishes indeed the harmonicity 
conditions. 

V. CONCLUSION 

In conclusion, this work considers the general post-Newtonian approximation formalism 
for isolated, slow-moving matter systems. Our motivation is to justify the practical 3.5PN 
calculations of radiation reaction achieved in paper I (Ref. Q). Our study follows from the 
previous investigation of paper II, in which a general scheme was proposed for the formal 
iteration of the post-Newtonian series in the source’s near zone up to any post-Newtonian 
order. There, the expression comprised: (i) a particular solution to the flat wave equation, 
given in terms of a (well-dehned) regularized solution, made with a hnite part FP, to the 
iterated Poisson equation for a non-compact source, and (ii) a homogeneous solution, dehned 
as a specihc multipolar anti-symmetric wave, associated with radiation reaction effects. 
Here, we present an alternative, complementary expression for the gravitational held in the 
source’s near-zone within the post-Newtonian approximation; notably, as in paper II, the 
expression incorporates fully the effects of the radiation reaction. It involves both particular 
and homogeneous solutions. In contrast to paper II, however, the particular solution in this 
instance represents a formal expansion of the retardations of the standard retarded integral, 
where each term has been suitably regularized using the FP procedure. The homogeneous 
second term, given in the form of a multipolar anti-symmetric wave, describes the non-linear 
contribution in the radiation reaction force of gravitational wave tails beginning at 4PN order 
~ 1/c®. The “ordinary” radiation reaction effects are then contained in the odd-parity part 
of the retardations of the particular, inhomogeneous solution. 

Consequently, as shown in Ref. |]|, the radiation reaction effects in the local dynamics 
of an isolated source up to and including the 3.5PN order ~ 1/c^ can be determined solely 
by considering the odd-parity part of the retardations in the standard retarded integral, 
following the appropriate regularization based on a FP prescription. 

Finally and complementary to the proof given in paper II, we verify that our post- 
Newtonian solution is divergenceless by demonstrating that the harmonic-coordinate con¬ 
dition is satished. Hence, the new expression is indeed a solution of the Einstein held 
equations. 
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Together with useful technical relationships given throughout this paper, the new form for 
the metric in the near-zone enables signihcant insights into the character of the gravitational 
radiation reaction and contribution of tails for isolated matter systems in the post-Newtonian 
expansion of general relativity. 
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